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Page  1,  line  3 from  bottom:  should  read  p . = ... 

n / jc 

2 

Page  1,  bottom  line;  should  read  E(z  .)  < “»  and... 

X f t 

Page  2,  line  3 from  bottom:  should  read  "where  = ••• 

Page  3,  line  8 from  bottom,  insert  comma  between  B and  B . 

Page  3,  line  4 from  bottom,  insert  comma  between  B and  B. . 

m 

Page  4,  line  11  from  bottom:  should  read  "converge  on  x S.," 

i=-m 

Page  4,  line  10  from  bottom,  add  "S.  ={F  ; |F  |<1}. 

1 ^ 

Page  5,  line  4 from  top;  should  read  E [ z „ . , a . ] 

q r 

Page  5,  Eq.(2.3)  should  read  Ynn  = I I -k- 

n=-p  k=l 

Page  6,  Eq.(3.1),  remove  parenthesis  from  exponent  e ^’^^^and 
place  below. 

Page  6,  lines  3 and  4;  e should  read  e . 

Page  7,  line  4,  beginning,  should  read  4>(B^,B^). 

Page  7,  line  6,  should  read  = (j)^^ 

Page  7,  line  2 from  bottom,  add  the  term  +y  ^ at  end  of  line. 

Z f di 

Page  8,  line  9 from  bottom,  should  read  Yqj^  = ‘J’2^-10*’' 

2 

Page  8,  line  3 from  bottom,  should  read  y _(00)  = o 

Z f di  cL 


Page  9,  line  4 from  top  should  read 


^0l”‘^2^-10"‘^1^00 

_2 


Page  9,  line  11  from  bottom:  should  read  Yn 

Page  9,  line  5 from  bottom.  Replace  "Taneja  et  al."  by  Voss  et  al. 

Add  the  following;  If  0^^  = @2  = 0,  the  ARMA  model  reduces  to  the 

1 dim  autoregressive  model  of  temporal  and  spatial  order  1.  The 
equations  for  the  autocovariance  then  agree  with  the  results 
obtained  by  Taneja  et  al. 

(over) 
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Page  2 


r 


Page  10,  line  5 from  bottom,  last  term  of  equation  on  that 

line  is  z . 

X,  t-2 

Page  11,  lines  3 & 4 from  bottom,  insert  "1"  at  end  of  each  line, 
Page  11,  bottom  line,  equations  should  read; 


0^  + ej  < 1 


®4  - S 


< 1 


e. 


64!  $ 1 


Page  12,  line  8 from  bottom,  last  term  in  equation  is  o^. 

Page  14,  line  4 from  top,  last  section  should  read 

{ (<l>4“04)  + 

Page  14,  line  5 from  top  should  read: 

Page  14,  line  6 from  top  should  read; 

,^^(-1-1)  - <VS2>'’a 

Page  14,  line  10  from  bottom  should  read; 

^01  “ *^1^00^*^2^-10'^*^3^-l-l'^‘*’4^0-l"®l®a"®3^‘^2"®2^°a”®4 '‘^I'^l^'^a 
Page  14,  line  9 from  bottom  should  read: 

Page  14,  line  8 from  bottom  should  read: 


11 


‘*’l^l0‘^‘^2^00‘*’‘^3^0-l'^‘^4^1-l"®2°a"®3  ^‘’*l"®l^®a 


Page  17,  lines  7,  6,  5,  4 from  bottom,  i should  be  replaced  by  <, 

Page  18,  lines  4,  5,  6,  7 from  top,  i should  be  replaced  by  <• 

Page  19,  line  7 from  top,  last  term  should  read  (000) 

Page  21,  line  9 from  bottom  should  read  p(g,f)  = p (g^^  ,g2 , f ) . . . 
Page  21,  line  8 from  bottom  should  read  p(gj^,g2,f) 

References,  second  page,  add;  Taneja,  V.A. , Aroian,  L.A.  (1977) . 
Time  series  in  m dimensions,  autoregressive  models. 
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1. 


INTRODUCTION. 


We  describe  a general  linear  stochastic  model  which 

supposes  a time  series  to  bo  generated  by  a linear  aggregation 

of  random  shocks  at  various  temporal  and  spatial  locations. 

Letting  x = (x,,x-,...,x  ),  a M dimensional  vector,  the 
12  m 

general  autoregressive-moving  average  model  (ARMA)  of 
M-dimensional  time  series  is: 

00  oo 

(1.1)  ^x,t  ~ n=-«)  k=l  l^^n,  k^x+n,  t-k  ®n , k^x+n , t-k]  ^ ^x,t, 

OO 

where  n = (n,  , n., , . . . , n ) and  7 denotes  the  M-dimension  sura 
12  m ^ 

n~“oo 

00  GO  oo  OO 

7 7 7 7 

over  ouch  of  the  M components  of  n (i.e.  n^=-«“'n  =-«^ 

1 2 m 

and  z ^ = z ._-E{z  ^),  the  deviation  from  the  mean. 

X^t  X/ti.  x^t 

The  a are  independent  random  shocks,  so  that  their 

X / t 

autocovariance  function  is:_ 


Y . = E ( a 

’n,k  X 


ff  n=o,k=o  i.e  n = (o,o,.... 


o,  otherwise 


and  their  autocorrelation  function  is: 


1,  n=o,k=o 


otherwise 


also  assume  that  z is  a weakly  stationary  process,  i.e. 


E(J„  . and  | | , | 


UUu 

iwrr.r;.- 
jBsn:,:-  ■ 


OtsL  AWi'J’  rZ'.ynrit'Liui 


OO 


OD  i 
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Tn  this  paper  we  explicitly  focus  our  attention  on  the 
special  case  of  (1.1)  in  which  only  a finite  number  of  the 
coefficients  are  non-zero,  that  is: 


q r vs 

*x,t  " n=-p  k=l  '^n,k^x+n,t-k  n=-u  k=l  ®n,k*x+n,t-k^^x,t 


Whether  or  not  the  coefficients  are  zero  it  is  easy  to 
represent  the  process  (1.2)  in  terms  of  shift  operators. 


'X/ 

% 

®t*x,t 

II 

N 

X 

rt 

1 

>-> 

B z 

= z f J 

X£  x,t 

x-6.' 

'V, 

x,t 

x+6  ^ , 

where  , • • • ' ‘im>  >'>'3 


1/  i = j 

(^otherwise 


Powers  of  these  operators  are  defined  in  the  usual  manner, 
for  example : 


2 '\,  'Xj  % 

B Z = B (B  z ^)  = Z 
X^  x,t  x^  x,t  x-26 


i,t 


-1 


In  addition,  we  note  that  B is  the  inverse  of  F i.e.  B = F 

x^  x^  , x^  x^ 


*1  -'i » 

Equation  (1.2)  can  be  written  in  terms  of  these  shift  operators: 
a r vs 


V'J-P. 

k^l*^n,k  ^x 

k 

®t  *!z,t  ■ 

lie. 

n=-u  k=l  n,k 

where  = 

(f"i  f''2 

. . . ,F”m) 

= (B“"i  b"'^2  . . 

X 

’‘l,  ^2, 

^m 

1,  ^"2, 

.-n 


Equation  (1.3)  can  be  rewritten  in  the  form: 
q.  r V 8 


(1.4)  (1-  ^ ^ (b  F^n^iz  - n-  ^ ^ ® 

n— p k=l  '*’n,k^x®t^*x,t  ~ n«-u  k-1  n,k^x\^®x,t 
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♦(B  ,B^)z  ^ 

= 0(B^,B.)a 

X t X , t 

X t 

<I>(B  B.) 

= -KB  ,B„ 

t 

Xi  x^ 

e(B  B) 

= 0(B  ,B 

*1  ‘ 

Xi  X2 

m 


X '“t' 

m 


We  may  consider  this  as  an  ARMA  model  of  temporal  orders 
r and  s,  and  of  spatial  orders  (p+q)  and  (u+v) . By  spatial 
order  p+q  we  mean  of  order  dimension  i for  each 

of  the  M dimensions  of  the  autoregressive  portion  of  the  model. 
Similarly  for  the  moving  average  portion  of  the  model,  the 
spatial  order  is  u+v  or  each  dimension  i.  We 

then  refer  to  this  as  a M-dim  ARMA  model  of  order  (r,s;  p,q;  u,v) 

The  ARMA  model  of  equation  (1.4)  may  be  considered  as  a 
M-dim  autoregressive  process: 

= ®x,f  ®x,t  * ‘"®x'®t’“x,t 

is  a M-dim  moving  average  process. 

This  model  may  also  be  considered  as  a M-dim  moving 
average  process: 

z . = 0(B  B )b  , where  b follows  the 

XtXft  X/C 

M-dim  autoregressive  process  defined  below: 


♦'“x'®t’'x,t  - 0<Bx-\>»'°x'®t>'’x,t  ■ 0<®x''t'“x,t 


4 


The  moving  average  terms  on  the  right  of  equation  (1.4) 
will  not  affect  the  conditions  for  stationarity  of  an  auto- 
regressive process  in  M-dimensions . Therefore,  these 
restrictions  on  the  parameters  of  the  autoregressive  model 
alone  will  apply  to  the  combined  ARMA  model.  The  condition 
for  stationarity  is  that  the  autocovariance  generating  function 


♦ (11^,0^)  must  converge  for 


[B^I  £ 1 and 


B < 1, 

X — 


B I 1 i = 1,2,...  ,M. 

^i 


we  mean 

.Similarly  the  conditions  for  invertibi lity  of  the  process 
that  apply  in  the  moving  average  model  also  apply  in  the  ARMA 
model.  The  condition  for  invertibility  is  that  n(B^,B^)  = 

0 must  converge  on  SqX.S_j^xS_2X.  . .xS_^,  where 

= {®x.‘  l®xj  - Sq  = {®i=  l^tl  1 4' 


2 . AUTOCOVARIANCE . 


The  autocovariance  of  the  mixed  ARMA  process  may  be  gotten 

% 

by  multiplying  equation  (1.3)  by  t-K'  '^^ere  I = (K,j^,2,2,. 

and  taking  expectations. 

Writing  equation  (1.3)  as: 


. , ) 
' m 


'V, 

z 


r — 

r 

<(>  . 

f"b'"z 

V s 

- 1 1 

0 ^ 

F^B^a  , 

n=-p 

k=l 

^n,k 

X t X , t 

n=-u  k=l 

n,k 

X t X , t 

q 

1 

n=-p 

r 

k=l 

*^n,k 

% 

z 1 

x+n, t-k 

V s 

- ^ 

n=-u  k=l 

0 . 
n,k 

^ x+n , t-k 

+ a 


Multiplying  by  Zx-ll,t-k' 

q r 
I J. 


V 


LI  1 L Z + 

x-i,t-ic^x,t  “ n=-p  k=l  '*’n,k  ^x+n,t-k^x-ll,t-K~n=-u  k=l  ‘^n,k^x+n,t-k^x-?,,t-ic‘'’ 


a ^ 

X,t  X-I!.,t-K 


Takinq  cxpoctaL ions ; 

q r 


V s 


''''11,1c  I ‘^n,k  ~ ^ ®n,k 

n--p  k=l  ' n=-u  k=i 


K-k)  + y _(^/<) 

Z f u 


■ '^'^x*n,t-k  Vt,t-K  is  the  croks  cxkrari^ 

'X/ 

between  z and  a and  is  defined  by:  y^  ~ x-£,t-  ,\,t 

'V' 

Since  z 4.  depends  only  on  shocks  up  to  time  t-x, 

X X>  ^ t ” K 

it  follows  that  y (P.,k)  = 0,  k > 0.  Since  the  a are 

Zd  X / u 

independent  random  shocks  it  also  follows  that  =0,  > 0 . 

Z a 

Consequently ; 

(2-2)  Y,,,  = I J *„,k  ''t+n,  k-k'  ^ 

rp-p  k=l 

**'i  ^ u^+1,  for  some  dimension  i,  i = 1,2,...,M. 

The  variance  of  this  process  is: 

<2-2)  Yoo  = I Y„,_k  - i ! e„,k  Y,,3(n,-k)  t oj 

n=-p  n=-u  k=l 

For  a given  ARMA  model,  this  system  of  equations  can  be 

solved  to  determine  the  parameters  <(>  , and  9 , in  terms  of 

n ,k  n,k 

the  autocovariances  y 

n,K 


3 . Spectrum: 

The  spectrum  of  a mixed  process  in  M-dimensions  follows 
from  tht  ,'i-dimensional  autoregressive  model  and  the  M-dimensional 
moving  average  model  in  a manner  similar  to  the  zero  dimensional 


case  of  Box  and  Jenkins. 


1 


(3.1)  p(g,t)  = 2<>^  I 0 (e  > | ^ 

It  ( o’^'Ce- 


0 £ q.  £ *5. 
0 < £ < \, 


where 


g = ,g2, . . -gj^)  and  i = and: 


0(e-2’i‘',<!-i2,a,  . „ ^^-i2„g^_^-i2i.g2 ^-i2.g____ ^-12™  £, 

♦ 3-121.0,  , t(e-i2''d,...,3-i2"g^_3-i2.£| 

4 • Partial  Autocorrelation: 

The  mixed  ARMA  process  may  be  written  in  the  form 
^x,t  ^ ^^x'^t^^^^x'^t^^x  t'  ® 

infinite  series  in  B ,B  ,...,B  ,B  as  in  the  moving  average 

Xi  X2  Xj^  t 

M-dimensional  model  and  therefore  the  partial  autocorrelation 
is  quite  similar  to  that  of  the  pure  moving  average  process. 

5.  ARMA  Model  (l.li  1,0;  1,0) 


We  consider  the  1-dim  mixed  process  of  temporal  order  1, 
and  spatial  order  1,  for  both  the  autoregressive  and  moving 
average  portions. 

From  equation  (1.3)  we  obtain  the  following  model: 


0 1 


0 1 


'5-11  ^X,t  ' ♦n,k  < \ ^x,t  - * "x.t 

*^-11  Vl,t-1  ^ *^01  ^x,t-l  ■ ®-ll  ®x-l,t-l  " ®01  \,t-l  \,t 
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From  equation  (1.4),  we  qet  the  followinq  form: 


For  this  model  then: 


For  convenience  let: 


4'0  ■ <>i 

iLil  ~ <(’2 


®0  " ®1 
0_11  = 02 


The  model  then  may  be  written  as: 

% '\/ 

z . = <b-z  , . , + 6,z  ^ ,-5_a  , . .-0,a  ^ , +a 

x,t  2 x-l,t-l  ^1  x,t-l  2 x-l,t-l  1 x,t-l  x,t 

From  the  corresponding  first  order  1-dimension  auto- 
regressive model,  this  ARMA  model  is  stationary  if  |4i2^  + 4>2l  ^ ^ 
and  l4ij^-((i2l  < 1 » or  equivalently  + |(|)2l  < 1*  From  the 

corresponding  1-dim  first  order  moving  average  model 
this  ARMA  model  is  invertible  if  19^^  + 02  1 1 and  |62”e2l  I' 

or  equivalently  |9^|  + I92I  < !• 

The  autocovariance  function  of  this  model  may  be  obtained 
from  equation  (2.1): 


8 


The  variance  for  this  first  model  is: 


(5.3)  Y 


00 


From  equation  (2.2)  , for  H 2 , or  < ^ 2 the  autocovariance 
depends  only  on  the  autoregressive  coefficients  4)^^  and 
We  therefore  obtain  the  following: 

''e.K  ■ ' V-i-iYVo-1 

The  remaining  autocovariance  terms  depend  on  both  '^l'*^2 

and  from  equation  (5.2)  are  given  as: 

(5.5)  - ♦jY.i  * »l''00-“2''za<-l°’  - ®l’'z,a'“'»^  ^z,a«>l' 

■-lO  ' <'2''0-J*Vl-rVza'‘>-l>-  “ll'z,a'^‘l>  * I'z.a'l'” 


’'ll  = '''2’'00+h’'l0-  Vza"’"’  • ‘’l’'z,a‘l'”  * ''za'l^> 

Since  this  ARMA  model  is  of  temporal  order  1 and  spatial  order  1 
for  both  the  autoregressive  and  moving  average  portions  of  the 
model,  the  only  non-zero  cross  covariance  terms  are: 

’'z,a"><»  - ”1 

’'z,a'“-»  - ‘♦r‘’l’  “a 

y a(-l-l>  = <*2-"2^  “a 


Substituting  these  cross-variance  terms  in  (5.3)  and  (5.5) 
yields  the  following: 


9 


^00  = °a^h^0-l-*-  ‘<’2^-l-l-  <'<'2"‘^2)^a 

- A .A  ..  n 2.  „ _ ^01"  2^-10''*’l^00 

^01  ^2^-10  ^1^00"^l‘’a'  1 ■ 2 


^10  “ ‘*’2^0-1  ‘^1^1-1  ‘ ‘*2  *^0 

ni  " <'2’'00*♦ln0-"2'’^■  »2  - liTliloOlllIlO 

a 

The  expressions  for  0 and  0^  may  be  substituted  in  the 
equations  for  Yqq  andYj^Q,  which  results  in  a pair  of  quadratic 
equations  in  4>  and  that  may  be  solved  for  ((>  and 

and  subsequently  for  0 and  0^  in  terms  of  the  autocovar iances 

y, 


Hrr- 


If  i))^  = 4>  ^ ~ 0,  the  AFIMA  model  reduces  to  the  1 dim 


moving  average  model  of  temporal  and  spatial  order  1.  The 
equations  for  the  autovariances  then  agree  with  the  results 
obtained  by  Taneja  et  al. 

The  spectrum  of  this  mixed  first  order  process  is  obtained 
from  equation  (3.1)  as: 

2o2|e(e“^2iTf^^-i2Tigj|2 


p(f,g)  = 


(e 


-i  2nf  -i2ng.  1 2 


,_2|.  « -i2Tif  -i2TTf  -i2Tiq|2 

!. ’ , 0 < f < S 

li  A -i27Tf  , -i27if  i2iiq|2  „ , 

|l-((ij,e  -ip^o  e-  ^1  0 1 g i ‘j 
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As  in  the  case  of  the  autocovariances,  if  6^^  = ©2  = 0, 
the  spectrum  reduces  to  that  of  the  1-dim  first  order 
autoregressive  model;  and  if  <p^  = 4)^  = 0,  the  spectrum  reduces 
to  that  of  the  1-dim  first  order  moving  average  model. 


6.  ARMA  Model  (2,2;  1,0;  1,0): 

The  next  model  considered  is  a 1 dim  mixed  process 
that  is  of  second  order  in  time  for  both  the  autoregressive 
and  the  moving  average  portions  and  of  first  order  in  the  space 
dimension  for  both  the  AR  and  MA  portions,  i.e.  an  ARMA  model 
of  the  form  (2,2;  1,0;  1,0). 


From  equation  (1.3)  we  obtain  the  following  form  for 
this  model: 

0 2 _ i.-v  0 2 


z 


. = y )■(()„  .-  y ye  KP^Bj^a  .+a  . 

;,t  r^-1  k=l  *^'^  ^ ^ n-lk-1  ^ ^ 

(6.1)  = ♦ 


,n„k. 


-l,l*x-l,t-l'^‘^-12,  ’'x-l,t-2'^‘‘’0l^x,t-l'^‘*’02  ^x,t-z 


~®-ll  ^x-l,t-l~®-12  ^x-l,t-2~®01^x,t-l~®02®x,t-2'''^x,t 


From  equation  (1.4)  we  get  the  following  form: 
^^"'‘’Ol^t"*^-!!  ®x®t”‘^-12  ^x®t"**'02®t^  ^x,t  = 

^^■®01®t"®-ll  ®x®t“°-12  ®x®t-®02®2)^x,t 
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For  convenience,  let: 


01 

^1 

4>. 

-11 

♦2 

4>. 

-12 

^3 

*^02 

'*’4 

01 

°1 

0_ 

■ 11 

®2 

0_ 

■12 

«3 

®02 

®4 

The  model  then  may  be  rewritten  as: 


% 


x,t  ‘^l^x,t-l’^‘^2^x-l,t-l‘^‘^3^x-l,t-2'^‘^4^x,t-2~®l®x-l,t-l 


~0-a  , . ,~6,,  a , o~0,a  . -+3 

2 x-l,t-l  3 x-l,t-2  4 x,t-2  x,t 


From  the  corresponding  1 dim  AR  model  of  temporal  order  2, 
and  spatial  order  1,  this  mixed  ARMA  model  is  stationary  if: 


+ 4'2  't’3  *l’4l  1 ~ <l>2  '<’3  “ '<'41  ^ ^ 

l4>3^  + 412  “ ^3  “ ^41  ^ 1 - <l'2  “ 4>3  + 04I  < 1 


From  the  corresponding  1 dim  MA  model  of  temporal  order  2, 
and  jpatial  order  1,  this  mixed  ARMA  model  is  invertible  if 

the  roots  of  = d"  ( 0i+®2®x^ ®t"  ^ ®4'^®3®x^  ®t  outside 

the  region  |B  | <l  and  |B  | < 1,  therefore: 
t ^ X . ^ 


02+  63+  0^  < 1 

01-  02- 

03+  0 

02+  O3+  O4  < 1 

O1+  ®2~ 

03+  0 

10  + ,0  1 <1 

I84-  63I  < 1 

^ 1^31^4! 

< 1 
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The  rtutocovariance  function  Cor  this  model  may  be  obtained 
from  equation  (2.1): 

" '*’-ll'^jl-l,K-l  ‘^-12^Jl-l,tc-2  ‘^Ol'^Jl,ic-l  ‘^02^J!„K-2''®-ll^za^‘^"^''^"^^ 

■°12^za^^'^''^-2)  - ■®02^za^^''"'2)  + 

= ^iY£,k-1  ^ ‘*’2Vi,k-1  ^ ^3Vl,<-2  ^ U'^Z,k-2 


The  variance  for  this  ARMA  model  of  temporal  order  2 and 
spatial  order  1,  is: 

(6.3)  Yqq  - + ^2Y-1-i  <>3Y-1-2  ^4V2  ' Vza^®"^^  ‘ Vza‘'^~^^ 

- 03V‘-1-2>  - Oiy^^O-2)  * yl 


From  equation  (2.2)  the  autocovariance  function  for  i ^ 2 
or  K > 3 only  depends  on  the  autoregressive  coefficients 

*^1'  *^2'  *^3'  *^4‘ 


We  then  have: 


“ '^-ll^ii,-l,K-l  * ^-l2^l-l,K-2  * ‘^01^£,k:-1  ‘^02^)1, k-2 

* £ iC“l  ^ ^2^£  — 1 iC”!  ^ *^3^£*“1/K“2  ^ ^4'^£/K'"2 


for  £ > 2 or  K > 3 


The  remaining  autocovariance  terms  depend  on  both  the 
autoregressive  and  the  moving  average  coefficients  and  are 
given  as  obtained  from  equation  (6.2): 

(6.5)  + ♦jY.j.i  ♦ ♦jY;,,!  - OjY^(0,0)  - SjVjal-l'Ol 

no  ■ nn-i  * Vo-i  * nva  * V1-2  ■ 

- 03Y„«.-2)  - 6,Y,a(l.-2)  * 


ni  “ ♦ii'io  * nno  * nn-i*  n^'i-i  ■ ni'za'i''”  - nna"’-"’ 

- Vza">'-“  - Vza<l--«  * na<l'“ 

I'-U  = V-10  * V-20  * ♦3l'-2-l  * V-1-1  ■ "ina'-^-'"  ' 

- «3-'za'-2'-»  - nna'-l'-”  * na'-^'l> 


^2  ’ ♦l''ol  * ♦2''-ll  * n^io  ♦ nno  - - nna'-'-'* 

- 633^ (-1.0)  - 03Y„(0,0)  . y^(0,2) 


112  ■ ♦ii'ii  * Voi  * nno  * nno  ■ nna'i’«  • nna'"-'' 

- nv^a"’-'”  - Vza'l'‘»  * Iza'^-^’ 


since  this  Al<i4A  model  is  of  temporal  order  2,  and  spatial 


order  1 for  both  the  autoregressive  and  moving  average  portions, 
the  only  nonzero  cross  covariance  terms  are: 


j(0  0) 

2 

= a 

a 

(4>4“64)  + 0^ 

^(0-1) 

Y,a(-l-2)  = 

'^1^V®2^  + (j)2(<(>j^-ej^)  + (<1)3-03)  a‘ 

^(-1-1) 

= (()>3^-e)o^ 

Substituting  these  cross  covariance  terms  in  equations 
(6.3)  auxi(6.5)  yields  the  following; 

^00  “ ‘^4''^0-2  " '’a“®2 ^‘^2~®2^ °a 

- 03{(()j^((|)2-02)+4)2(<('i-0i)  + (<(>3-03))o3  - e4{  (^4-04)+^! }o^+o^ 

^01  • '^iW‘<’2Y-10^<'3Y-1-1^4Vi®1  a " ®3^'^2"®2^''a  -®4^h-®l^®a 

YlO  * ♦l^l-l"’Vo-l^Vo-2‘'Vl-2  ■ 02^V®l^V®3^<V®4^^'^l^'<’r®l^^^a 

^11  * ‘^lYl0^♦2^00^‘^3^-l♦^^l-^®2'’^®3^'^^®l^‘’a 

Y-ir  nY-10-^‘^2^-20^^3Y-2-l^‘^4^-l-r®4^‘^2"®2^^a 
^02  = nY01^^2^-11^^3^-10^^4^00-°4‘^a 
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^lYll-'*<’2Yoi^^3YoO^‘^4^10‘®3°a 


Through  an  iterative  procedure  this  set  of  AR  and  MA 
coefficients  may  bo  solved  in  terms  of  the  autocovariances. 

If  = <(>2  = 't>3  “ 'l’4  “ 0»  this  ARMA  model  reduces  to  the  1-dim 
moving  average  model  of  temporal  order  2 and  spatial  order  1. 
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The  above  equations  for  the  autocovariance  function  agree  with 
the  results  obtained  by  Voss  et  al.  If  the  ~ ~ 

are  zero,  the  model  reduces  to  the  1-dim  autoregressive  model 
of  temporal  order  2 and  spatial  order  1.  These  autocovariances 
agree  with  the  results  obtained  by  Taneja  et  al. 

The  spectrum  of  this  mixed  process  of  temporal  order  2 and 
spatial  order  1 is  obtained  from  equation  (3.1)  as: 

/ c\  o 2 I -i27rg  -i27rf.  i2 

p(g,f)  =2^1  0(e  ^,e  ) 1 

I -i2TTg  -i2Tif)  i2 

I <Ke  ^,e  I 

- 0^2  I , - .-i2'irg  -i2Trf  ^-i2iTg  -i2TTf  ^ -i2irf  -i4TTf  i 2 

= 2a^  I 1-02©  ^e  -O^e  ^e  -G^^e  I 

I , . -i2irf  . -i27rf  -i2Trg  . -i4iTf  -i2'ng  . -i4Trf  1 2 3 

I l-<t>2e  e e \ , | 

0 < f < »5 

0 ^ g £ *5 

As  in  the  case  of  the  autocovariance  function  if 
= <p2  - 4>2  ~ 'f’4  ~ spectrum  of  the  mixed  process  | 

reduces  to  that  of  the  1-dim  moving  average  model  of  jj 

temporal  order  2 and  spatial  order  1.  If  0^^  = ^2  ~ ^2  ~ '^4  ~ ^ ' I 

the  spectrum  reduces  to  that  of  the  1-dim  autoregressive  model  1 

of  temporal  order  2 and  spatial  order  1.  j 


r 
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7 . AUMA  Model  (1,1;  1,0,  1,0;  1,0,  1,0) : 

The  next  model  considered  is  a 2 dimensional  mixed 

model  of  first  order  in  time  and  each  of  the  2 space  dimensions, 

for  both  the  autoregressive  and  the  moving  average  portions, 

i.e.  and  ARMA  model  of  the  form  (1,1;  1,0,  1,0;  1.0, 1,0). 

X dim  y dim  x dim  y dim 
AR  AR  MA  MA 

From  equation  (1.3)  we  obtain  the  following  form  for 
the  model: 

- y y y 0 .F^IF  2 a . + a . 

- 1 -_i  v-1  n,,n-,)c  X y t x,y,t  x,y,t 


1 2 


% 


^ 'J’-l-ll  ^x-l,y-l,t-l  *^-101  ^x-l,y,t-l  *^0-11  ^x,y-l,t-l 

* *^001  ^x,y,t-l  ^-1-11  *x-l,y-l,t-l  °101  ^x-l,y,t-l 
*^0-ll  ®x,y-l,t-l  ^001  ^x,y,t-l  ^ ^x,y,t 


From  equation  (1.4)  we  get  the  following  form: 

'*’001  \ “ '^-101  ®x®t  " '^0-11  Vt  ' *^-1-11  ®x®y®t^^x,y,t 
“ ®001  \ " ®-lbl  W “ ®0-ll  Vt  " *^-1-11  ®x®y\^®x,y,t 


I 

I 


r 

I 


; i 

\ < 


! 

I 

t 


For  this  model  then: 
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1 

t 

I 


♦'“xVt’  - 1 
‘"'“xjVt'-  ^ 


'♦ooi'* 


(6  -9 

' 001 


-101 

-101 


B -<b 
X 


B -6 
X 


0-ll®y*<t>-l-ll 

0-ll®y"®-l-ll 


B B )B 
X y 


B B )B 
X y 


t 

t 


For 

convenience , 

let: 

*^001 

= h 

'<’-101  = '<'2 

'<’0-11  = '<’3 

'<’-1-11  = 

%01 

«-l0l  = «2 

Vll  ^ ^3 

0-1-11  = 6, 

The  model  then  may  be  rewritten  as: 


^x,y,t  “"^l^Xjy.t-l  * ‘^2*x-l,y,  t-1  ^ '^3*x,y-l,  t-1  ^ *^4^x-l,y-l,  t-1 


®l^x,y,t-l  ®2^x-l,y,t-l  ®3*x,y-l,t-l  ®4^x-l,y-l, t-1  ^ ^x,y,t 


From  the  corresponding  2 dim  autoregressive  model  of 
temporal  order  1,  and  spatial  order  1 in  both  x and  y directions, 
this  mixed  model  is  stationary  if: 

I <t>i  + <<>2  *^3  ‘**4  I ~ ^ 

I (j)j^  - (>2  + <l>3  - <l>4  I 1 

Ul  + 4'2  “ <<>3  ~ '^’4  I i 1 

I (j»j^  - ^2  ~ *^3  *^4  I ^ 

From  the  corresponding  2 dim  moving  average  model  of 

temporal  order  1,  and  spatial  order  1 in  x and  y,  this  mixed 

= 0~^(I!  ,B  , B ) converges 

A y t 


! 


i 


i 


i 

t 


j 


model  is  invertible  if  ir(B  ,B  ,B^) 

X y t 
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1 X S_  j X ^_2'  wJ>c*ro  ^ l|  and 

■i  ■ "'®x'By=t>=  [i-<VW«3VWy>®t]''- 


this  mixed  model  is  invertible  if: 


'®1  ■"  ®2  ®4'  - ^ 

I®1  - O2  ■ ®3  ^ ^4!  1 ^ 

|0^  + 02  - 83  - O4I  < 1 

I03  - O2  + 83  - ej  < 1 


The  autocovariance  function  for  this  model  may  be  obtained 
from  equation  (2.1): 


(7.2)  Y 


<,3,  *,2''^*  '’’-l-ll^)l3-l,£2-l»K-l'*’'^101^«,3-l,£2''«"l^'^0"l^^^l'J*'2"^''^~l 


lYd 


"®0-ll^za^*'l'*'2“^''^^^"®00l''^za'*'l'*'2''^'^^'^^za^^l'^2''^^ 

“"l^za  ^ *■!' *'2' ~^2''^za  ^ *"2' ’’*^3^za  ^^1' *'2~^' 

"®4^za  ^2“^' 


2a'''l"'2' 


The  variance  for  this  ARMA  2 dim  model  of  first  order  in 
time  and  in  each  space  dimension  is: 

*■'•  ’>  ^000  “ ♦l''oo-l  *V-10-1  * ♦3V1-I  * V-1-1-1  ■ “ll'za 
- Vza'-^"»  - Sl-za ">-!-«  - 


The  remaining  autocovariance  terms  depend  on  both  the  auto- 
regressive and  the  moving  average  coefficients  and  are  given  as; 


^>01  ~ *1^00  * *2^-100  * *3^0-10 

^4^za 

'^101  “ ^1^100  * ^2*000  * ^3^1*10  *^4^0-10  " 

- “2^4 WO'’*  - “sl-za'^-^'  - 04Vza'°-“>  * 

'^011  ■ *1^010  * ‘'’2^-110  '*'  *3^000  ^4^-100  ” 

-'>2^za<-l“'  -aj'.alOOO)  - 0^3 ^^(-100)  *3^(011) 

''ill  ' *l''U0  * *2''oiO  * '*3'' 100  '’'4''000  ■ “llza*^'''’* 

- 023„(010)  - 033,^(100)  - 0,3,a(000)  3 3^,(m) 

’100  " ♦l''l0-l  * ♦2''00-l  * ♦l^l-lO  * ♦4’'0-l-l  ■ ®l''za'“-“ 

- »2''za<“-»  - «3^za'l-l'»  ' Vza"’'^'”  ''za'*'’”' 

''010  ■ ♦l’'01-l  * ♦2’'-U-l  * *3''00-l  * V-lO-l  ■ ‘>l'>za"’^‘^’ 

- '’2''za<-^‘-“  - '’3''za"’'’-^>  ' ^''za''^-^’  * ''za''’^* 

''no  ■ ♦l''ll-l  * ♦2''01-1  * bl'10-1  * ♦4''00-i  * “Hza'll-” 

-0^3^^ (01-1)  -033^^(10-1)  -033^(00-1)  *3^(110) 

"l-lO  ■ *l''l-l-l  *■  *2'' 0-1-1  * "'3'' 1-2-1  * ♦4'0-2-l  ■ Ollza'^'^’'’ 

-6^3^(0-1-1)  -O33 ^^(1-2-1)  -0,3 ,,(0-2-1)  3, ^,(1-10) 


r 


20 


Since  this  ARMA  model  is  of  temporal  and  spatial  order  1, 
for  both  the  moving  average  and  the  autoregressive  portions, 
the  only  nonzero  cross  covariance  terms  are; 

Y^{000)  = ol  Y,^(00-1)  = ^,^(-10-1)  = [V02^^a 


Y,^(0-1-1)  = (,>.^-,.^10,^ 


Y^^  (-1-1-1)  = [^0410: 


Substituting  these  cross  covariance  terms  in  equations 
(7.3)  and  (7.5)  yields  the  following: 

^000  = ‘^1^00-1  ^ '^^2^-10-!  ^ ^3Vl-l  ^ ^4Y-1-1-1  " 

^001  ^ h^ooo  '^2^-100  ‘‘’3^0-10  ^ ‘^4^-l-10  " ^I'^a 

^101  ^ '^1^100  '*’2^000  ‘^3^1-10  ‘^4^0-10  “ 

^011  ^ '^1^010  *^2^-110  *^3^000  '^4^-100  “ '^3'^a 

^111  " '*’1^110  * VolO  ‘^3^100  '•'4^000  ” °4‘^a 

^100  " '•'l^lO-l  ^ ‘*’2^00-1  * '‘’3^1-10  Vo- 1-1  " ^2^'^r^^l^°a~^4^‘^3‘®3^°a 

^010  = '^1^01-1  ‘^2^-ll-l  ^ VoO-1  ^ V-10-1  ■ - ®4<V®2^°a 

''^110  Vll-1  * ‘^2^01-1  VlO-l  '^4^00-1  “ ”4^'^r°l^°a 

2 

"^1-10  ” ^^l"^!- 1— 1 ^^2''^0— 1— 1 *^3'^1— 2— 1 ^4^0— 2—1  **  ^^2  ^a 


Through  an  iterative  procedure,  the  coefficients  for  the 
AR  and  MA  portions  may  be  solved  in  the  terms  of  the  autoco- 
variances. If  <t>2^  = <('2  ~ <l>3  = <(>4  = 0,  this  ARMA  model  reduces 
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to  the  corresponding  2 dim  moving  average  model  of  temporal 
and  spatial  order  1.  The  above  results  for  the  autocovariance 
function  agree  with  the  results  of  Voss  et  al.  If  the 

= Oj  = ''3  = = 0,  the  model  reduces  to  the  corresponding 

2 dim  autoregressive  model  of  first  order  in  time  and  each  of 
the  space  dimensions.  The  autocovariance  then  agrees  with  the 
results  of  Taneja  et  al. 

The  spectrum  of  this  mixed  process  of  dim  2,  and  first 
order  in  time  and  each  space  dimension  is  obtained  from 
equation  (3.1)  as: 
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